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Abstract 

Ground state properties and excited states of ladder-type paraphenylene 
oligomers are calculated applying semiempirical methods for up to eleven 
phenylene rings. The results are in qualitative agreement with experimental 
data. A new scheme to interpret the excited states is developed which re- 
veals the excitonic nature of the excited states. The electron-hole pair of the 
Si-state has a mean distance of approximately 4 A. 
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I. INTRODUCTION 



The investigation of 7r-conjugated polymers is in many ways affected by the structural 
disorder in these systems. In contrast, the ladder-type poly-p-phenylenes (LPPP) |J offer the 
opportunity to study large, rod-like chains of planarised phenylene units. As a consequence, 
the 7r-system might spread out over an entire polymer and a vibronic resolution of the 
Si <— So transition is discernible ||. In order to deduce some characteristics of the polymeric 
films P], like the effective conjugation length, several oligomers have been synthesized in the 
past to study the low- lying electronic excited states of the polymer . 

Photoconduction in LPPP films |J has been measured as a function of the energy of 
the exciting light, too. A typical small plateau of the photocurrent occurs between the 
absorption edge and its steep rise at higher energies and extends in this case over 1.6 eV. 
This behavior of the photocurrent which does not follow directly the absorption features is 
sometimes called "abatic". 

One possible explanation for this effect rests on the interpretation of the electronic ex- 
cited states for the individual molecules. Excited states of 7r-conjugated molecules are usu- 
ally described as Coulomb bound electron-hole pairs. This physical picture originates from 
solid-state physics of (organic) semi-conductors. Therefore, these molecular states are often 
referred to as excitons, although they have to be clearly distinguished from the extended 
band states in a crystal. 

A reasonable estimate of the exciton binding energy in conjugated polymers has been 
determined, e.g., by scanning tunneling spectroscopy measurements [0] which have lead 
to a value of about 3.5 eV. Excited states with a smaller value, and larger electron-hole 
distance, respectively, should be more susceptible to the separation via an external electric 
field. Following this idea, the conjecture has been brought forward that such a state is 
responsible for the steep rise of the photocurrent in poly-phenylene-vinylene (PPV) ||. 
Later on, another explanation has followed based on the excess photon energy which is 
converted to the vibrational heat bath M. The latter proposal is now widely accepted. 

In order to test these concepts for excited states of 7r-conjugated systems, several 
oligomers of the LPPP type with up to eleven phenylene rings are investigated in this 
work. The study of oligomers instead of an (infinite) polymer follows the above mentioned 
approach and allows the direct comparison with experiment. The main difference to the 
experiments in condensed phases is the restriction to single chains in the vacuum. 

As the experimentally used molecules are computationally too demanding one has to 
replace the large aliphatic substituents attached to LPPP by hydrogen (see Figure [I] and 
Table |). This should have only negligible effects on the optical properties, which are gov- 
erned by the frontier orbitals of 7r-symmetry. These aliphatic substituents are only necessary 
for the better solubility of the polymer, or to prohibit the formation of aggregates in the 
film (R3 = methyl in Figure |TJ). 

Since the systems studied here reach the size of the effective conjugation length proposed 
for LPPP (about 14 phenylene rings [0,0), ab-initio or density functional methods are not 
applicable, and one has to assent to less expensive semiempirical methods (AMI, INDO/S; 
see below). Thus, the wave functions of the corresponding ground states are the INDO (in- 
termediate neglect of differential overlap) Slater determinants |$ NDO ) [[10]. For the excited 
states 



$s ) (see equation [16]) , the INDO/S expansion is used in the spin-singlet sector. The 
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excited states with dominant oscillator strengths will be addressed as Si for the first excited 
state, S m for the intermediate excited state and Sf for the high energy, "Frenkel-type" state. 
The electronic ground state will be denoted as S = |$o NDO ). 

The article is organized as follows. In Sect. |I|, the semiempirical approach is briefly 
described. In Sect. fTJ, the results for the geometric structure of the oligomers and their 
spectra are presented. The main part of this article, Sect. [IV], focuses on the development of 
a general interpretation scheme for excited-state wave functions. Its application to INDO/S 
wave functions leads in a straightforward way to the interpretation of the excited states as 
bound electron-hole pairs. A short conclusion closes the presentation. 



II. METHODS 

Although it is not feasible to calculate the higher oligomers by first-principle meth- 
ods, the oligomer with two phenylene rings (n = 0.5) has been calculated at MP2/6-31G* 
level HTI| , |l2| (M0ller-Plesset Pertubation Theory of second order). The results are used as a 



reference for the semiempirical methods. 



Following a procedure of Bredas et al. [13|, one has to choose a semiempirical method 



which yields reliable geometric structures. In the present case the Austin semiempiri- 
cal parametrization (AMI) [jl4[ gives an almost exact agreement with the MP2 results 
(cf. J2 I A r c-c| for the bond lengths in Table Hj). This method will therefore be used to calcu- 



late the geometric structure of the ground states. Note, however, that the PM3 method (15 
yields better results for zero-point vibrational energies (ZPE). 

The AMI and MP2 calculations have been carried out on a IBM/SP2 computer using 
the GAUSSIAN94 (Rev. D4) |16| program package. All minima are characterized by diago- 
nalization of the Hessian Matrix, whereby the zero-point energies (ZPE) have been obtained, 
too. The latter will be given unsealed throughout the paper. 



In the next step, excited singlet states are calculated using Zerner's INDO/S method fl7 
based on the minimum AMI structures from the electronic ground states. Thus it is clear 
that geometric relaxation effects in the excited state are neglected. The active CI space 
consists of the 22 highest occupied and the 22 lowest unoccupied molecular orbitals. This is 
the biggest active space possible within the used program package and it contains orbitals 
of a symmetry for some oligomers as well. As expected, the dominant configuration state 
functions in the wave function describe n-n* excitations. The calculation of the spectra 
have been accomplished on a Pentiumlll-PC using the program package CAChe 3.1 JTS 



In order to get a more realistic view of the calculated line spectra, Gaussian peaks are 
least-square fitted to the INDO/S oscillator strengths. This procedure masks transitions 
with moderate oscillator strengths which are close in energy to a dominant transition. This 
is especially the case for the high-energy region of the spectrum, in which strong Sf <— So 
transitions cover others. The hidden states are not visible in the optical spectrum, and 
consist of a number of important configurations sensitive to the size of the active CI space 
(Figure ^|). The remaining part of the spectra comprises only features with one transition, 
except in two cases, where two almost identical states are close together. In these cases only 
one state will be discussed (the Si <— So transition for the n = 0.5 oligomer and the S m <— So 
transition for the trimer; see discussion below). 
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III. RESULTS FOR GEOMETRY AND OPTICAL SPECTRA 



A. Geometry 

All structures have been optimized without any symmetry constraints. Therefore, the 
geometries only fit to certain point groups within crystallographic accuracy (± 0.001 A, 
±2°). The resulting point groups are equivalent to ones suggested by the simple Lewis- 
type structure, see Figure [TJ. Oligomers with an even number of phenylene rings adopt Civ 
symmetry. Those with an odd number adopt C^h symmetry. The C-C bond lengths range 
between about 1.380 A and about 1.500 A. Figure |] shows the MP2 result for the n = 0.5 
oligomer as an example. 

As the smallest C-C bond distance in a methylene bridge is 1.503 A long, they are 
assigned to single bonds. As a consequence, the hydrogen atoms of the methylene bridge 
do not participate in the conjugated system, although this would be allowed. This supports 
the validity of neglecting the aliphatic substituents. 

B. Optical spectra 

In Figure § the (vibronically unresolved) calculated spectra of LPPP oligomers are plot- 
ted. A comparison with measured fluorescence spectra of different oligomers [|5] shows qual- 
itative agreement, i.e., they show a broad Si <— So transition which is shifting to lower 
energies with increasing system size. As one would expect for the Si state, the HOMO 
LUMO excitation is the dominant one. HOMO (LUMO) refers to the highest occupied 
(lowest unoccupied) molecular orbital in the INDO ground state. 

A second interesting feature can be seen in the calculations for the oligomers with five 
and more rings. It will be referred to as the S m <— S transition, because it is energetically 
in the middle between the Si state and the high energy region. In the S m state the dom- 
inant determinants are built by the substitution of the HOMO by the LUMO+2 and the 
substitution of the HOMO— 1 by the LUMO+1. The position of the S m <— So transition 
also shifts to lower energies with increasing oligomer length. This is not discernible in the 
experimental spectra because of the small corresponding oscillator strengths. 

Finally, both sets of spectra show a steep rise at around 5.40 eV which is dominated by 
a transition at approximately 5.85 eV for every oligomer. This transition will be called the 
Sf So transition, as the Sf state shows a high degree of localization (i.e., a "Frenkel" 
state), which will become clear later in the discussion. The Sp state is determined by several 
configuration state functions, where low-lying occupied molecular orbitals are exchanged by 
orbitals with high energy. 

In Figure [| the energy of all three optically important transitions is drawn as a function 
of the reciprocal number of phenylene rings in the molecule. The transitions of the Si 
and the S m states show a strong linear dependence, whereas the energy of the Sf <— S 
transitions hardly changes with system size, in agreement with the experimental observation. 
Despite this qualitative agreement, a quantitative comparison with the experimental values 
for the effective conjugation length [[| is not possible, since the theoretical curve leads to an 
unphysical negative value. The main reason for this error is the neglect of the polarization 
energy in the calculations which stabilizes the excited states in the condensed phases in the 
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experiments. Moreover, the linear extrapolations for the transitions to the S m and the Si 
states suggest a crossing of the two energies which must not be taken for granted at this 
stage. 

From the electronic dipole transition moments M for the Si <— So transitions, one obtains 
an estimate for the radiative fluorescence lifetimes r for every oligomer, using 



M = e( $ 



CIS 



$ INDcj 



3/ic 3 



T 



3I M |2 



(la) 
(lb) 



Here, e is the elementary charge, fi is the position vector of the i-th particle of charge ^e, 
h is Planck's constant, c is the speed of light, and hu is the energy of the transition, all 
quantities measured in cgs units. The lifetimes r show a linear dependence on the reciprocal 
number of phenylene rings, see Figure |^. Here, the value of several hundred picoseconds 
for the large oligomers is in quantitative agreement with the experimental value of about 
300 ps for the polymer, measured in an organic matrix [[TjJ. Under the assumption of a 
linear relation between r and the inverse of the oligomer length, one finds from the data 
in Table |TV| that a value of r po i ymer = 300 ps corresponds to an effective conjugation length 
of about 20 phenylene rings, in reasonable agreement with the value of 14 rings estimated 
from other optical experiments 



IV. INTERPRETATION 



A. General considerations 



For an interpretation the wave functions for the ground state and the excited states need 
to be related to experimentally observable quantities. The optical absorption at frequency u 
is proportional to the real part of the optical conductivity, as given by the Kubo formula 



Re [a(u > 0)] 



Im [ Xjj {uj > 0)] 



(2) 



where Xjj( u ) is the current-current correlation function, 

1 



Xjj(uJ > 0) 



E 



E — H + u + irj" 

l<*. lil^o)! 2 



$ 



u - (E s - E ) + ir\ ' 



(3) 



Here, H is the Hamilton operator of the system, are its eigenstates with energies E s 
(s = 0, 1, 2, . . .), h = 1, and 77 = + is positive infinitesimal. Therefore, the real part of the 
optical conductivity may be written as 



Re [a(u > 0)] = - £ |($, \j\ $ )| 2 5 (u - (E s - E )) 



(4) 
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The current operator is defined by 



Pn 



(5) 



where Cl a , C n ,a are creation and annihilation operators for electrons with spin a =f, { 
in the molecular orbitals n , and p min is the matrix element between the corresponding 
one-particle states. 

Equation (f|) is readily interpreted. The absorption of a photon with energy u induces an 
excitation between the ground state |$o) and the excited states |$ s ) with energy u) = E s —E . 
The amplitude for this absorption process 

Am = <$ s |5| $ ) (6) 

2 

determines the oscillator strength, / 0jS cx |A 0iS | . 

For a further analysis of the amplitudes A 0tS , the current operator is expressed in terms 
of the field operators, 

*i(3=£#(2)<5k> (7a) 

n 

*<r(»)=Z)0n(5)C'n,<T, (7b) 
n 

which create/annihilate an electron with spin a at x. From the completeness of the molecular 
orbitals (p n (x) one readily obtains 

j=Y,l dxdyj{x, y)i>t(x)it a {y) (8) 

with 

j(?,y) = ^Vm,n<l>n{x)(j)* m {y) . (9) 
n,m 

Therefore, the amplitudes can be cast into the form 

A o,s = 53/ d %edx h j(x e ,x h )A% s (x e ,x h ) , (10) 

where the particle-hole amplitudes 

A° s (x e ,x h ) = ($ s ^(£^04) $\ (H) 



are given by the overlap between the excited state |$ s ) and the ground state |$o) with an 
electron at x e and a hole at (if x e ^ Xh)- 

It is thus seen that the quantities Aq s (x e , Xh) allow to address the question in how far a 
given excited state |<J> S ) may be viewed as an electron-hole excitation of the ground state, 
and they directly enter the Kubo formula for the optical conductivity. Note that the analysis 
of (|TT|) does not require the full spatial dependence of the many-particle wave functions |$ s ) 



or |$ ) as in |8|j21|] or the density-density correlation function of the ground state as in p2| . 
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The matrix elements for the current operator, j(x e ,Xh) in @, do not change much over 
atomic distances. Therefore, it is usually sufficient to introduce the coarse-grained densities 
for the electron-hole content of |$ s ) with respect to |$o)> 



PoA^j) = J2j d %edx h A% s (xe, %h) ©(^ - ri)Q(x e - fj) 



(12) 



where 



0(x - fj) 



1 , if x - n e Vi 
, else 



(13) 



is the step function for the (atomic) volume Vi around atom i. P 0tS (i,j) is the overlap density 
between the excited |$ s ) and the ground state |$o) with an electron-hole pair around the 
nuclei at i and j, respectively 

Equation (|TT| ) is exact and applies to all quantum-mechanical systems. It is currently 
under investigation for the numerical analysis of quantum lattice systems with the density 
matrix renormalisation group (DMRG) method [23|. In the next subsection, the present 



case of INDO/S wave functions will be studied in more detail. 



B. Application to INDO/S 

As described earlier, the ground state is approximated as an INDO Slater determinant 



l*o> 



INDO 




INDO 

cr,0 



ni $ T°)=n n ci> ac ™>> cu*) 

a a l<a<n 

— J dx% . . . dx n Det [4> a {xi)] ^t(^i) • • • ^i(^n) (vacuum) , (14b) 



where the molecular orbitals 4> a (xi) (1 < a < n) are expressed as a linear combination of 
spatial atomic orbitals (AOs) Xb{%) 



AOs 

U^) = da ( b )Xb(xi 

b 



(15) 



The atomic orbitals are centered at certain nuclei such that b = (fb, (3b) comprises the 
orbital type = s,p, d . . .) and its position r&. 

The spin singlet excited-state wave functions are described as linear combinations of 
singly-excited configuration state functions 



$cis 



$ CIS 

s 



INDO 




(16a) 
(16b) 



where the indices a and r refer to the active space of occupied and virtual molecular orbitals, 
respectively. 
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With these approximations for |$o) and |$s), it is obvious that the coordinates of the 
added electron enter (11) only through the virtual MOs (j) r whereas the hole is described by 
the eliminated, occupied MOs 4> a , 



INDO\ 
/ 



(17) 



To make further progress, one has to address the coarse-grained densities Po :S (hj) m (|12|), 



^indo/s 

0,s 1 



P, 



INDO/S / . 



0,s 



hJj 



dxpdxi 



^Car(t)* r {x e )(t)a{x h ) 



9(xh ~ ri)Q(x e - fj) 



(18) 



E E<rC„ v K(r J> /3 J ))V(r J! /3 J )(d a '(f i ,A))Wi ! A^ • 



In the last step the INDO approximation |]T0[ and the orthogonality of the atomic orbitals 
on the same atom were used, i.e., 



(19) 



Recall that d r (fj,Pj) denote AO coefficients, and c Qjr are the CI coefficients of the CIS wave 
function. In this case, one may verify that J2i,j Ps,o(h j) — 1j if one takes into account the 
INDO approximation and the normalization of |$P S ) in flllf) , 



< ISi^CIS^ __ J 



J2c* r (j)*(x e )(j)a(x h ) 



nuclei 



dx P dxh 



E ^ 



INDO/S/ • 



s,0 



(20) 



h3 



Therefore, P^ D °^ S (i, j) may be viewed as a discrete probability function for finding a hole 
around the nucleus i and an electron around the nucleus j. This simplification is only valid 
in the INDO/S approximation. 

A check can be made, whether the values of P™ DO//S (i, j) carry significant weight corre- 
sponding to AOs of a symmetry. It turns out that this weight is vanishingly small, and one 
can sum all contributions of the AOs of an atom without losing information. 

One has to verify that the detected electron-hole correlations for an excited state really 
stem from the electron-hole interaction rather than from a coincidence of the motion of two 
independent particles confined to a small molecule. Therefore, the quantity 



lNDO^ 







(21) 



has been investigated according to the scheme (|TT[) . Here, |$h^l) is the INDO wave function 
where the HOMO is substituted by the LUMO. This corresponds to the pure "band-like" 
excitation, i.e., in a CIS state representation this configuration state function has the coef- 
ficient c ar = 1 for a = HOMO, r = LUMO, and all other coefficients are set to zero. 
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Given the P s {hj) an expectation value (E) of the distance of the electron-hole pair 
can be calculated with the help of the bond lengths r(i,j) = \f\ — fj\ of the molecule, 

(E) s ,o = Y,pT° /S ^J)^,j). (22) 

The standard deviation <t Sj o is equally accessible, 

Ko) 2 = £(r(i,j) - (E) S , ) 2 P™ DO/S (i,j) . (23) 

y 

In small molecules, the confinement of the oligomers will determine the distance of the 
electron-hole pair. In larger molecules, if the interaction between hole and electron is weak, 
(E) S fi and cr Sj o will increase with system size since the particles move essentially indepen- 
dently through the molecule. On the other hand, the interaction between electron and hole 
may keep them together at a fixed distance even though the size of the molecule increases. 
Such a bound electron- hole pair may lead to a constant value of (E) 8j q and a Sj o for every 
system size. Note, however, (E) Sj o and ct S; o only contain an information about the overall 
extension of an excitation but may fail to describe the localization of an electron-hole pair 
onto segments of the molecule; an example of this situation is given below. 

In order to get a more pictorial way of the electron-hole pair distribution, one may 
concentrate on a quasi one-dimensional chain of carbon atoms; Figure |] shows how those 
chains are chosen. When P St0 (i,j) is plotted for this chain a bound electron-hole pair will 
show large values along the diagonal of the plot and vanishingly small values in the off- 
diagonal regions. Unbound pairs will lead to the opposite situation. 

The two-dimensional distribution P Sj a(i,j) can be further smoothed into 

/-r+Ar/2 nuclei 

P.,o(r) = dr'Y, P sAhj) 8 (r' - \n - r,\) , (24) 

Jr-Ar/2 i ■ 

which is solely a function of the electron-hole distance. The choice of Ar = 1.8 A gives 
smooth curves as a function of r. P 8 ,o(r) gives the most concise description of the electron- 
hole excitation. 



C. Results 

The -P s I o D °^ S (^i) matrix has been calculated for the states, whose transitions are of 
significant oscillator strength in the spectra: The Si, the S m , and the Sf state. 

First, the state Si is discussed. As seen from Figure 0, no bound state is formed for 
only two phenylene rings. In fact, the molecular confinement is dominant up to a chain 
length of four phenylene rings. For five phenylene rings, see Figure [8], and larger systems, 
a bound electron-hole pair is discernible. The shape of Ps u o(i,j) in the plot along the one- 
dimensional chain does not change for systems larger then five phenylene rings. The Si state 
clearly corresponds to a strongly bound electron-hole pair. 

Figure |9| shows the average electron-hole separation and the corresponding standard 
deviation for the Si state as a function of system size. They saturate for more than five 
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phenylene rings, (E)s 1 ,o ~ 4 A. This is in reasonable agreement with the experimental value 
of 7 A obtained from electro-absorption measurements [[|. Some results of electron energy- 
loss spectroscopy, however, seem to hint at a totally different behavior of the excited states: 
they are supposed always to increase with increasing length of the oligomer P^j . At present, 
the reason for this discrepancy is not clear. 

The saturation behavior for (E)$ 1 q and the apparently large saturation value for the 
standard deviation o"g l ~ 3 A are readily understood from the "smoothed" probability 
distribution Psi,o( r ), see (pif). As seen in Figure |10|, Psi,o( r — > 0) decreases with increasing 
system size, the whole distribution broadens, and develops a maximum around r m = 4 A. For 
oligomers with five and more phenylene rings the distribution does not change significantly. 
Since at the same time the system grows, the values of (E)s 1>0 and ctsi.q saturate. The 
relatively large values of &Si,o are due to the location of the maximum of the distribution at 
a finite value of r m . Finally, Figure 11 shows the probability distribution P™q°^ S (z, j) for 
the largest oligomer with eleven phenylene rings (n — 5). 

Next, the state S m is addressed. Bredas and coworkers have located a state in the 
calculated spectrum of poly-p-phenylene-vinylene (PPV), which they assign to a charge 
transfer state. In this context this state is equivalent to an unbound electron-hole pair. The 
optical transition to this charge transfer state is in the same region as the steep increase 
of the photocurrent in the respective polymeric film. As a result, they conclude that the 
population of this state is responsible for this "abatic" behavior of the photocurrent. The 
S m state obtained in the present calculation lies in the same region as the charge-transfer 
state in PPV, that is between the Si and the high-energy states. The question is, whether 
the S m state is a charge transfer state or not. 

As seen in Figure [T2], the S m state shows almost constant values of (P) Smi0 and o"s m ,o- 
A rise in these values occurs for the n = 5 oligomer which hints at unbound states for 
larger systems. This idea, however, is not supported by a more detailed look at the graphs 
for Ps m ,o( r ) as a function of oligomer length, see Figure [L3|, and the results for Ps m ,o(z,j) 
for the largest oligomer, see Figure [14|. The curves are qualitatively the same as for the 
corresponding Si states. For example, compare Figures |TT| and |I4|: P Sl (i,j) and P Smi0 (z, j) 
are essentially zero in the off-diagonal regions. Similarly, the curves Ps m ,o( r ) as a function of 
oligomer length resemble those of Ps 1; o( r )> compare Figures [10] and [13]. The only difference 
to the Si state is the localization of the excitation to three parts of the molecule in S m . 
Because of this, one cannot interpret this state as an unbound elctron-hole pair or exciton, 
and consequently no explanation for the rise of the photocurrent of LPPP at about 4.0 eV 
can be given at this level of theory. This is in line with the assumption that local heating 
due to excess energy is the reason for the behavior of the photocurrent || and not the nature 
of the optically accessible states. 

Lastly, for the Sf state, the mean electron-hole separation does not saturate but grows 
with system size, see Fig. |15[ From such an analysis one might conclude a delocalized state. 
On the contrary, for the Sf state the degree of localization is higher than in the other two 
excited states, and the electron-hole pair is actually restricted to every single phenylene 
ring. As can be seen from Figures [16] and [T7|, the apparent lack of convergence in (P)s F ,o 
and <7s F ,o relates to the fact that isolated elctron-hole pairs on chain segments can be linearly 
superimposed, i.e., an extended state can be formed which, nevertheless, does not contribute 
to the conductivity. As for the states Si and S m , there is no spreading of Ps F ,o{i,j) into 
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off-diagonal regions with increasing system size, and the smoothed distribution function 
Ps F ,o( r ) displays the same trends as before. 

So far, every state investigated can be regarded as a bound electron-hole pair with 
different degrees of localization. They are molecular analogues of the excitons in the physics 
of semiconductors. In order to see the differences to a pure "band excitation" , the analysis 



is repeated for |$h->l), see (pip . As expected, and confirmed in Figures [L8[ 0, and |2(J, the 
motion of electron and hole in the molecule is uncorrelated. In contrast to the excitonic 
cases before, the almost linear increase of the mean electron-hole separation and its variance 
as a function of system size in Figure [T8| is accompanied by a broadening and flattening of 
the smoothed distribution function -PH^L,o( r ) ; see Figure |19[ As seen from Figure ^ there 
is considerable weight in the off-diagonal region in the probability distribution Pa^L,o(i, j), 
which actually looks like a half sphere. Hence, the excitonic behavior of the states Si, S m , 
and Sf is genuine, and not just a coincidence in the uncorrelated motion of an electron and 
a hole in a restricted geometry. 



V. CONCLUSION 

The ground and singlet excited states of various oligomers of LPPP have been described 
with semiempirical methods. A qualitative agreement was achieved with the experimental 
absorption spectra, especially for the measurements of the fluorescence lifetime, and the 
mean electron-hole distance. A new analysis of the excited states has been given. As in 
optical absorption experiments, the overlap matrix elements between excited states and 
the ground state with an electron-hole pair are studied as a function of their respective 
positions on the oligomers. Excited states with high oscillator strengths are found to be 
bound electron-hole pairs. Therefore, no explanation of the abatic onset of the photocurrent 
in LPPP films can be given at this level of a microscopic theory. 
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FIGURES 

FIG. 1. The LPPP oligomers; the substituents used by various experimental groups are sum- 
marized in Table |. In this work, all substituents are replaced by hydrogen. 

FIG. 2. Optimized geometric structure of the n = 0.5 oligomere of LPPP on MP2/6-31G* level 
of theory. 



FIG. 3. Calculated absorption spectra as a function of the oligomer length. 
FIG. 4. Calculated energetic position of the optical transitions depending on the system size, 



compare Table IV. The experimental values for the Si <— So transition are taken from 



FIG. 5. Calculated fluorescence lifetimes as a function of the system size, compare Table IV 



FIG. 6. For plotting P s $(i,j) along a quasi one-dimensional chain, only those carbon atoms 
are into account which are marked with bold lines. 

FIG. 7. Psi,o(^j) = P-ij for the Si state of the oligomer with two phenylene rings (n = 0.5). 
The state is determined by finite-size effects of the molecule, i.e., there is no bound electron-hole 
pair. 

FIG. 8. Ps lt Q{hj) = Pij for the Si state of the oligomer with five phenylene rings (n = 2). A 
bound Frenkel exciton is discernible. 

FIG. 9. Expectation value (E)g 1 fi and standard deviation osx.o f° r the state Si in units of A 
as a function of the oligomer length. 

FIG. 10. Smoothed distribution function Psi,o( r ) = Pij as a function of the separation r for 
the state Si for various oligomer lengths. 

FIG. 11. Probability distribution -Psi,o(«>j) = Pij along a quasi one-dimensional path for the 
oligomers with eleven (n = 5) phenylene rings for the state Si. 

FIG. 12. Expectation value (£ ; }s m ,o and standard deviation 0"s m ,o for the state S m in units of A 
as a function of the oligomer length. 

FIG. 13. Smoothed distribution function -Ps m ,o( r ) = Pij as a function of the separation r for 
the state S m for various oligomer lengths. 
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FIG. 14. Probability distribution Ps m ,o(*,j) = Pij along a quasi one-dimensional path for the 
oligomers with eleven (n = 5) phenylene rings for the state S m . 

FIG. 15. Expectation value (E)s F ,o and standard deviation <rs p ,o for the state Sf in units of A 
as a function of the oligomer length. 

FIG. 16. Smoothed distribution function -Ps F ,o( r ) = Pij as a function of the separation r for 
the state Sf for various oligomer lengths. 

FIG. 17. Probability distribution Ps F ,o(i,j) = Pij along a quasi one-dimensional path for the 
oligomers with eleven (n = 5) phenylene rings for the state Sf- 

FIG. 18. Expectation value (£ : )h->l,o and standard deviation cjh^l,o f° r the state |<3?h^l) i n 
units of A as a function of the oligomer length. 

FIG. 19. Smoothed distribution function PH^L,o( r ) = Pij as a function of the separation r for 
the state |<3?h^l) for various oligomer lengths. 

FIG. 20. Probability distribution Ph— >L,o(hj) = Pij along a quasi one-dimensional path for the 
oligomers with with (a) two (n = 0.5) and (b) eleven (n = 5) phenylene rings for the state I^h^l)- 
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TABLES 



TABLE I. Aliphatic substituents as defined in figure |]. 





1, § 1 


i], [|] this work 


Ri 
R 2 
R 3 


n-Hexyl n-Hexyl H 
p(n-Decyl)-phenyl p(t-Butyl)-phenyl H 
Methyl H H 


Ri 


H if terminus, Ri otherwise 



TABLE II. Comparison between MP2, AMI, and PM3. 



Method 


E / a.u. a 


ZPE/(kcal • moL 1 ) b 


point group 


E|Ar C -c| c 


MP2/6-31G* 


-499.8413498 


117.80 


Civ 




AMI 


0.0862897 


121.55 


Cl v 


0.000 


PM3 


0.0778001 


117.65 


Civ 


0.029 



a E = electronic energy (MP2), heat of formation (AMI, PM3) 
b ZPE= zero point energy 

c J2 \^ r C-c\ = Sum of the differences of the C-C bond lengths with respect to the optimized 
structure at MP2 level 



TABLE III. Ground state energies of the LPPP oligomers at AMI level. 



N a 


E/a.u. b 


ZPE/(kcal • moL 1 ) c 


E + ZPE /a.u. 


point group 


2 


0.0862897 


121.55 


0.279986 


Civ 


3 


0.1373472 


178.64 


0.422035 


Cih 


4 


0.1883913 


235.71 


0.564020 


Civ 


5 


0.2394355 


292.76 


0.705985 


Cih 


7 


0.3415265 


406.85 


0.989898 


Cih 


9 


0.4436183 


520.95 


1.273804 


Cih 


11 


0.5457102 


635.03 


1.557695 


Cih 


a N 
b E 


= number of phenylene 
= heat of formation 


rings 







c ZPE= zero point energy 
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TABLE IV. Calculated energetic position of the transitions from the ground state in cm -1 to 
the Si, S m , and Sp state, depending on the system size, and fluorescence lifetimes depending on 
the system size. The experimental values for the Si <— So transition are taken from pi. 



1/re a 


Experiment [5 


Si 


S m 


Sp 


t/10~ 12 s b 


0.5U0 




o A T A r o 




47 (03.1 


1 /o9 


0.3 


29500 


31214.5 




47629.4 


1683 


0.250 




29289.2 




47145.2 


1115 


0.200 


25700 


28055.1 


3861.4 


46995.0 


920 


0.143 


23800 


26841.0 


34345.4 


46851.2 


693 


0.1 




26187.2 


31613.2 


46758.9 


560 


0.09 




25795.0 


30272.1 


46686.4 


470 


0.083 


22100 











a n = number of monomer units 

b r = fluorescence lifetime following equation ([!]) 



TABLE V. Values of the expectation value {E) s ,o and the standard orientation a s ,o for the 
optically detectable states, depending on the size of the oligomers. 













N a 








state 


values b 


2 


3 


4 


5 


7 


9 


11 


|$h^l) 




2.93 


3.86 


4.76 


5.64 


9.01 


9.07 


10.79 




O"H^L,0 


1.80 


2.49 


3.16 


3.82 


5.77 


6.38 


7.65 


Si 


(^S!,0 


2.34 


2.96 


3.37 


3.58 


3.79 


3.87 


3.90 




CT Si,0 


1.13 


1.91 


2.35 


2.59 


2.82 


2.91 


2.96 


Sm 


(^)s m ,o 








3.06 


3.05 


3.20 


3.52 




°"S m ,0 








2.18 


2.19 


2.21 


2.59 


S F 


( E )s F ,o 


1.99 


2.12 


2.74 


2.64 


3.90 


3.59 


3.80 




O"S P ,0 


1.33 


1.49 


2.09 


2.03 


3.56 


3.11 


3.39 



a N = number of phenylene rings 
b [(E) )=[a]= 1 A 



17 



FIG. 1. Jorg Rissler Phys. Rev. B 
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FIG. 2. Jorg Rissler Phys. Rev. B 




FIG. 3. Jorg Rissler Phys. Rev. B 
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FIG. 4. Jorg Rissler Phys. Rev. B 





FIG. 6. Jorg Rissler Phys. Rev. B 
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FIG. 7. Jorg Rissler Phys. Rev. B 
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FIG. 12. Jorg Rissler Phys. Rev. B 
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FIG. 13. Jorg Rissler Phys. Rev. B 
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FIG. 14. Jorg Rissler Phys. Rev. B 




FIG. 15. Jorg Rissler Phys. Rev. B 
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FIG. 17. Jorg Rissler Phys. Rev. B 




FIG. 18. J6rg Rissler Phys. Rev. B 
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FIG. 19. Jorg Rissler Phys. 
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